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We propose a method to control the optical transparency of a Bose-Einstein condensate with 
working energy levels of the A-type. The reported effects are essentially nonlinear and are considered 
in the framework of an exactly solvable model describing the interaction of light with a A-type 
medium. We show how the complicated nonlinear interplay between fast and slow solitons in the 
A-type medium points to a possibility to create optical gates as well as to a possibility to store 
optical information. 
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I. INTRODUCTION. 

Recent progress in experimental techniques for the co- 
herent control of light-matter interaction opens many op- 
portunities for interesting practical applications. The ex- 
periments are carried out on various types of materials 
such as cold sodium atoms p], Q , rubidium atom vapors 
HHHH, solids 0,0, photonic crystals @. These ex- 
periments are based on the control over the absorption 
properties of the medium and study slow-light and su- 
perluminal light effects. The control can be realized in 
the regime of electromagnetically induced transparency 
(EIT), by the coherent population oscillations or other 
induced transparency techniques. The use of each dif- 
ferent materials brings in specific advantages important 
for the practical realization of the effects. For instance, 
the cold atoms have negligible doppler broadening and 
small collision rates, which increases ground-state coher- 
ence time. The experiments on rubidium vapors are car- 
ried at room temperatures and this does not require ap- 
plication of complicated cooling methods. The solids are 
obviously one of the strongest candidates for realization 
of long-living optical memory. Photonic crystals provide 
a broad range of ways to guide and manipulate the slow 
light. 

In this paper we study the interaction of light with 
a gaseous active medium whose working energy levels 
are well approximated by the A-scheme. Our theoreti- 
cal model is a very close prototype for a gas of sodium 
atoms, whose interaction with the light is approximated 
by the structure of levels of the A-type. The structure of 
levels is given in Fig. ^ where two hyperfine sub-levels 
of sodium state 3 2 Si/ 2 with F = \.F = 2 are associ- 
ated with |2) and |1) states, correspondingly [jj. The 
excited state |3) corresponds to the hyperfine sub- level 



of the term 3 2 P 3 / 2 with F = 2. We consider the case 
when the atoms are cooled down to microkelvin temper- 
ature in order to suppress the Doppler shift and increase 
the coherence life-time for the ground levels. The atomic 
coherence life-time in sodium atoms at temperature of 
0.9/xK is of the order of 0.9 ms 0. Typically, in the ex- 
periments the pulses have length of microseconds, which 
is much shorter than the coherence life-time and longer 
than the optical relaxation time of 16.3ns. The gas cell 
is illuminated by two circularly polarized optical beams 
co-propagating in the z-direction. One beam, denoted as 
channel a, is a c^-polarized field, and the other, denoted 
as b, is a cr + -polarized field. The corresponding fields are 
presented within the slow varying amplitude and phase 
approximation (SVEPA) as 



E = e a £ a e 



(k a z—uj a t) 



6b £bC 



i(k b z-ui h t) 



+ C.C. 



(1) 



Here, k a ^ are the wave numbers, while the vectors e a , e b 
describe polarizations of the fields. It is convenient to 
introduce two corresponding Rabi frequencies: 



_ 2p a £a „ 2fi b £ b 
S La — Z , "6 = Z 



(2) 



Here are dipole moments of quantum transitions in 
the channels a and b. 

Within the SVEPA the Hamiltonian describing the in- 
teraction of a three-level atom with the fields assumes 
the form 

H\ = Hq + Hi, fl =wia|l)<l|+w|3><3| ) (3) 
Hi = -\ {VL a e li - k - z ~ u ^\i){\\ + n b e t( - kbZ - UJ >' t '>\3)(2\) + h.c. 

Here H = 1. The description can be further simplified in 
the interaction representation for the density matrix p, 
namely, 
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In the new variables the Liouville equation takes the form 



|1) = |3S 1( ,,F = 2,M, =+1) 



\2) = \3S V2 ,F = 1,M T = - 1) 



FIG. 1: The A-scheme for working energy levels of sodium 
atoms. The parameters of the scheme are the following: 
lo 12 /(2tt) = 1772MHz, w/(2tt) = 5.1 • 10 14 Hz (A = 589nm), 
and A is the variable detuning from the resonance. 



Then the Liouville equation for the transformed operator 
p assumes the form 



where 



id t p= [(H -H E + Hj),p] , 



fl> = -£(fi«|3)(l|+n 6 |3)(2|) + fc.& 

o o n* a 
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n a n b o 



(4) 



The dynamics of the fields is described by the Maxwell 
equations 

(d 2 - c 2 d 2 z )Q a e^ z -^ = -^fl? (/5 31 e 4 ( fe ^-^*)) , 
{d 2 - c 2 dl)Sl h e l{ - kbZ -^ = ~^d 2 {p32e l{kbZ - Ubt) ) , 

where v a = {n A \p a \ 2 uj a ) / e , v b = (n A \p b \ 2 u b ) / n A is 
the density of atoms, and Cq is the vacuum susceptibility. 
For many experimental situations it is typical that the 
coupling constants v a ,b are almost the same. Therefore 
we assume that v a = v b = i/q. Within the SVEPA the 
wave equations are reduced to the first order PDEs: 



dtCla = Wo/031, d ( Q b = iv p32 



(5) 



Here C, — z / 'c,t = t — z / 1 c. 

Equations Eqs. 10 can be rewritten in a matrix form 

as 



d ( Hj=i^lD,p], 



(6) 



where 




drP 



—D 

2 



Hi,p 



(7) 



As was indicated above, the system of equations 
Eqs. ©,0 is exactly solvable in the framework of the 
inverse scattering (IS) method [Icj . 

This means that the system of equations Eqs. ©>0 
constitutes a compatibility condition for a certain linear 
system, namely 



9 C * = V(X) * 



» "OP iT; 
2 A— A w ' 



(8) 
(9) 



Here, A e C is the spectral parameter. The comparison 
against leads to the zero-curvature condition 
U c (A) - K(A) + [C/(A), V(X)] = 0, which holds identically 
with respect to the linearly independent terms in A. It is 
straightforward to check that the resulting conditions co- 
incide with the nonlinear equations Eqs. ©,10. At this 
point it is worth discussing the initial and boundary con- 
ditions underlying the physical problem in question. We 
consider a semi-infinite C > active medium with a pulse 
of light incident at the point £ = (initial condition). 
This means that the evolution is considered with respect 
to the space variable Cj while the boundary conditions 
should be specified with respect to the variable r. In 
our case we use as the asymptotic boundary conditions 
the asymptotic values of the density matrix as |r| — > oo. 
To solve the nonlinear dynamics as described by equa- 
tions Eqs . (JHJ) , (0> , the IS method considers the scattering 
problem for the linear system Eq.JSJ), while the auxiliary 
linear system Eq. 10 describes the evolution of the scat- 
tering data. The purpose of this work is, in particular, 
to study an essentially nonlinear interplay of the fields in 
the both channels. This goal leads to consideration for 
equation Eq. 10 the scattering problem of finite density 

type (cf. 01 an d references therein), i.e. H, a ,b —> ^1 as 
|t| — > oo. 

The scattering problem for equation Eq.0, reformu- 
lated as a matrix Riemann-Hilbert problem, is then posed 
on the two-sheet Riemann surface for the spectral param- 
eter A. The purpose of this work is to discuss certain ex- 
perimentally relevant solutions of the A-system, while the 
technicalities of the IS method are very well reflected in 
the existing literature anyway. For an account of results 
for the A-system accessible thr oug h the IS method see, 
for example, references [ill IT2I Il3l IT ^| . We only wish to 
indicate that in our analysis we use elements of the ap- 
proach based on the Riemann-Hilbert matrix problem, 
along with algebraic techniques of the theory of solitons. 
We use Darboux-Backlund transformations, in the spirit 

of [niii mm, up to certain modifications, however 

(cf. Appendix). The paper is organized as follows. In 
the next section we introduce the notions of slow and fast 
solitons and describe a nonlinear mechanism of formation 
of the transparency gate for the slow soliton. Section III 
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describes the transparency gate for the fast soliton. The 
section VI is devoted to conclusions. The details of the 
Darboux-Backlund transformation for the A-system are 
given in the Appendix. 

II. THE TRANSPARENCY GATE FOR THE 
SLOW SOLITON 

In this section we introduce a concept of slow and fast 
solitons in the A-medium and explain how the nonlinear 
interplay between the solitons leads to a possibility to 
control transparency of the medium. We discuss first the 
mechanism of transparency control for the slow soliton. 
We explain how the fast soliton propagating in the a 
channel hops to the b channel where the slow soliton is 
propagating. The fast soliton then "knocks down" the 
slow soliton, thus stopping the propagation of the latter, 
and then disappears itself due to the strong relaxation in 
the system. 

As was indicated above, in this work we consider exact 
solutions of the Maxwell-Bloch system Eqs.JJJiO exist- 
ing on a non- vanishing finite background. In our con- 
siderations the background field plays the same role as 
the controlling field in the conventional linear theory of 
EIT. The background field enters the exact solutions as a 



parameter in a substantially nonlinear fashion. Because 
of reasons of experimental relevance we specify the back- 
ground field as 

f^ 0) = cos(r?)n, ftj; 0) = sin(7 ? )f2, fl = n Q e lkC . (10) 

Therefore S7q specifies the total intensity of background 
fields existing in the both channels, r\ defines the relative 
intensities, while k <C k a ,b is introduced in order to take 
into account small spatial variations of the phase. The 
intensity of the background field Qq is an experimentally 
adjustable parameter. We show in the paper that this 
parameter provides the control over the transparency of 
the optical gates and determines the speed of the slow 
soliton in the system. 

For simplicity, we first assume k = 0. The atoms are 
initially unexcited and are prepared in the so-called dark- 
state 

|Vd)=cos(j7)|2)-sm(r7)|l) (11) 

Hereafter we assume r\ = 0, and therefore the dark state 
is simply |2). Taking this state of the atomic subsystem 
as the initial state we find the following solutions describ- 
ing the formation of the transparency gate (cf appendix 
for details): 



— — 



+2a2a3£o) 



.„, f "0*0 



fl h = -2i- 



/2^a ie ^^^(A-Ui^^V^-^y^ + ^/^^+^e-iV^-^y^,^/^^) 



>+a%e 



(12) 



The solution Ea. H12|) describes the nonlinear interaction 
of the slow and fast solitons. As is explained in the ap- 
pendix, the solution is parameterized by the constants 
defining the position and phase of the solitons. Without 
loss of generality we assume a\ and CI3 to be real con- 
stants, while a 2 = 1. We show below that a\ determines 
the position of the slow soliton whereas 03 determines the 
position of the fast signal. In practice these constants are 
defined by the initial condition, which specifies the actual 
pulse of light entering the medium at the point C = 0. 
To understand the structure of the slow soliton one can 
put 03 = 0. This choice corresponds to taking the fast 
soliton to —00 in the variable r. Indeed, this specifica- 
tion removes from the overall solution Ea. l|12l) the fast 
pulse component and thus singles out the slow soliton 
part. The slow soliton solution assumes then the follow- 



ing form: 



Cl a = £!otanh(0 s ) (13) 
cosh(0 s ) 



where 



c 



£0^0 



2(A 2 + £ g) 2 



- eo-Weg-ng +ln|ai|, (15) 



is the phase of the slow soliton. From the expression 

e 2 

above and in the simplifying approximation > > 1 the 
group velocity of the slow soliton can be easily derived: 

^(A 2 + £ 2) 



The pure state of the atomic subsystem corresponding 
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to the slow soliton solution Eg. l|14[) reads 



7^ + ^> |2> 



6 2Vah^ v / £0 - % /^! 



(16) 



Notice that the population of the upper level |3) is pro- 
portional to the intensity of the background field. The 
speed of the slow soliton is also proportional to Slg. This 
means that the slower the soliton, the smaller the pop- 
ulation of the level |3) and, therefore, the dynamics of 
the nonlinear system as a whole is less affected by the 
relaxation process. 

To understand the structure of the fast soliton one can 
choose a\ — 0. We then arrive at an expression describing 
a signal moving on the constant background with the 
speed of light (fast soliton): 



n i 



£0 



cosh(0/) 

where the phase of the fast soliton is 



O 6 = 0, 



(17) 



<j> f 



ill 



ln\a 3 \ 



We wish to emphasize once again that the atomic sub- 
system is prepared in the dark state |2), which is not 
affected by a signal propagating in the channel a. This 
explains why the fast soliton propagates with the veloc- 
ity of light. The fast signal propagating in channel a can 
be observed as a localized peak in the field intensity I a . 
Below we analyze how the fast soliton collides with the 
slow soliton, hops into channel b 7 and then slightly slows 
down. 

Figure[21illustrates the propagation and collision of the 
slow and fast solitons according to equation Ea. l(T2|l . The 
figure for I a shows the intensities of the signals in chan- 
nel a. The slow soliton corresponds to a groove in the 
background field £lo- It is clearly seen that after the act 
of collision the slow soliton ceases propagating in channel 
a, while some trace of the fast soliton still can be noticed 
in that channel. The figure for I a is complemented by the 
figure for the intensity lb of the field in channel b. We 
see that before the collision only the slow soliton exists in 
channel b, while after the collision the slow soliton disap- 
pears and a fast intensive signal appears, whose velocity 
is slightly below the speed of light. The process described 
above can be summarized as if the fast soliton knocks 
down the slow soliton. The notion of a transparency gate 
requires the existence of two distinctly different regimes, 
which are transparent (open gate), and opaque (closed 
gate) . In the absence of the fast soliton the gate is open 
for the slow soliton. When the fast soliton is present the 
slow soliton is destroyed, while the fast intensive signal 
created after the collision in channel b is attenuated due 
to strong relaxation in the atomic subsystem. The gate 
thus closes up in the course of the dynamics due to the 





s it 




FIG. 2: Knocking down the slow soliton. The two upper plots 
correspond to the dynamics of the intensities of the fields Q a 
and Sib. The two lower plots show the populations of the levels 
|1) and 1 3). The background field Qq = 1. Other parameters 
of plotted solutions are always chosen to be the same through 
the paper, namely ci = C2 = cz — 1, Ao = 2i, vq — 3, A = 0. 



relaxation process. To further explain this process we 
provide the Fig. [2] plots for populations of the levels |1) 
and 1 3). 

Notice that before the collision the population of the 
upper atom level |3) is negligible and is approximately 
given by the formula for the slow soliton solution Eq. i|16|) 
(see the lower right plot of P3). The populations of the 
lower levels |1, 2) are determined by the slow soliton (see 
the lower left plot of PI). Indeed, the fast signal existing 
in channel a does not interact with the atoms because 
at the onset of the dynamics their state coincides with 
the dark state |2). Figure |2 shows that after the collision 
the atoms of the active medium are highly excited and 
therefore the level |3) is strongly populated. This leads 
to the fast attenuation of the speedy intensive signal in 
channel b due to the relaxation. The optical gate closes 
up. 

To this point we have described a mechanism of con- 
trolling the transparency of the medium for a particular 
type of slowly moving signals. We now discuss a possi- 
bility to read information stored in the atomic subsys- 
tem. Let us assume that the background field vanishes, 
i.e. Hq = 0. As was explained above, the speed of the 
slow soliton then vanishes as well. However, the infor- 
mation about polarization of the slow signal is stored in 
the atomic subsystem. This effect can be interpreted in 
terms of the concept of a polariton, which is a collective 
excitation of the overall atom-field system. The notion 
of a polariton for the A-system has been used before. In 
the linear case the dark-state polariton was discussed in 
|19| . In the strongly nonlinear regime, which is the case 
for the present work, a similar interpretation is possible. 
Indeed, the field component of a slow soliton solution 



5 



can be interpreted as the light contribution into the slow 
polariton. When the controlling field f2 vanishes this 
contribution vanishes as well, along with the speed of 
the polariton. The latter then contains only excitations 
in the atomic subsystem. The general solution Eq. (|12|) is 
then reduced to the form: 



— 4icic 3 e exp sor 



2c 2 c 3 coshO s0 )+c? exp[2e O T- i 0^ 



C3 " 



(18) 



where 



3s0 



, 2 (A 2 +e a ) n ~ m ( c 2/c3) is the phase of the 
slow soliton for the vanishing background f2o- The form 
of the fields resembles a superposition of fast and slow 
solitons in Eq. (|12[1 , with the vanishing velocity of the slow 
soliton. It is impossible to single out the contribution of 
the slow pulse by choosing C3 = 0, because the amplitude 
of the slow signal approaches zero as the background field 
vanishes (see Ea.(|14p). The atomic state describing the 
stored information reads 



A+ie 
A— ie 



|2> 



(19) 



2^ 



-(*< 



ci 



1 2) + 1 3) 



For ci = the fields vanish, while the atomic state re- 
duces to a form corresponding to a stopped polariton 
described by Ea. (|16|l with Qo = 0. In other words, when 
the slow soliton is completely stopped its information is 
stored in the spin polarization of the atoms. As long 
as the upper state |3) is not populated, the state of the 
atomic subsystem is not sensitive to the destructive in- 
fluence of the optical relaxation processes. 

The conventional way Q to read the information 
stored in the atoms is to increase the intensity of the 
background field. Our method of reading the informa- 
tion is different. We propose to send the fast soliton into 
the space domain in the active medium, where the in- 
formation is stored. The polarization in the domain is 
then flipped by the fast signal. This is how the reading 
of information is realized. This way of reading optical in- 
formation is advantageous because it involves fast easily 
detectable processes. Figure illustrates the mechanism 
of the reading. Notice that the act of reading, based on 
the polarization flipping induced by the fast signal, can 
be realized in a very short time scale compared to typical 
relaxation times. 



III. THE TRANSPARENCY GATE FOR THE 
FAST SIGNAL 

In this section we explain how to create the trans- 
parency gate for the fast signal rather than for the slow 
soliton as in the preceding section. Instead of decreasing 
the background intensity we now propose to increase it 




FIG. 3: Reading the optical information by the fast soliton. 
The two upper plots illustrate the dynamics of the fields Q a 
and Qt- The two lower plots show the populations of the 
levels |1) and |3). The standing peak on the plot for Pi corre- 
sponds to the stored information in the form of the localized 
polarization. The rapidly moving localized excitation of the 
atoms given on the plot for P3 represents the act of reading. 
The background field Qq — 0. 




FIG. 4: Knocking down the fast signal. The small localized 
polarization flip is the slow soliton part. The intensive peak 
is the fast signal. The background field Qq = 1, Ao = i. 



to reach the point where f^o = So- As we show below (cf 
Fig.0]as well), the fast signal can be stopped by the slow 
soliton. In other words, the control of the transparency 
for the fast signal is realized by tuning the background 
intensity to the parameter eq entering the slow soliton 
solution. This effect is related to the appearance of the 
so-called exulton type solutions reported before , [2(j 
for the SIT and NLS models. 

For the constant background field with k — and for 
= £0 the solution has a simple and transparent form, 
namely 



-2e-*»i ((c 2 +c 3 tQ ) 2 -3cI) 



q n 4»e 1 e 2(n + A " ) (c 2 +c 3 (l+Tiio)) 

6 " C f e ^i+2 e -^i((c 2 + C3 mo) 2 +c§)' 



(20) 
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and 



A + iSlg 

A-if2 



|2) 



^ C2—C3+C3TQ0 
C2+C3+C3TQ0 



1) 



C2+C 3 +C 3 T(l 



|2) 



(21) 



-2(A-i!2 ) 



13} 



Here, the phase cf> s i corresponds to 4> s with Eq = Qq. 

To explain the mechanism underlying the functioning 
of an optical gate for the fast signal we analyze the solu- 
tion Ea. H2Q(l in detail. In the case C3 = we obtain the 



conventional slow soliton solution. If c\ — we arrive at 
a rational algebraic expression, which is a solution of the 
exulton type. In the case k = and c\ = the pulse, 
whose algebraic form is a rational expression (exulton), 
moves with the speed of light without decay. 



It can be readily seen that if k is not zero the exul- 
ton part of the solution vanishes in the course of the 
dynamics. For simplicity, and to demonstrate the typical 
behavior of this kind, we choose c\ = 0, c 2 = 0. Then, 



a, 



fi I 1 - 2e ifeC 



(iA(l + rOp) + n (l - ik( + rn ))(iA(l - rfip) + » (-l + ik( + tQ )) 
A 2 + (l + (fcC~ Ar) 2 )r>2+T 2 ^ 



n b = 0. 



(22) 



r 



The collision of the fast signal (exulton) with the slow 
soliton results in the disappearance of the former in chan- 
nel a. This means that the slow soliton controls the trans- 
parency of the medium for the signal moving with the 
speed of light. The reported mechanism is a realization 
of an optical gate for the fast signal. 



IV. 



CONCLUSIONS AND DISCUSSION 



In this paper we discussed a concept of the trans- 
parency gate for the fast and slow solitons in a A-type 
medium. We explained how the fast soliton can knock 
down the slow soliton and close the gate for the latter. 
The opposite process of closing the gate for the fast sig- 
nal, given by rational algebraic expression (exulton), is 
also described. We also described the process of read- 
ing optical information, written into the active medium 
by the slow soliton. It is worth discussing here a pos- 
sibility to actually create in the A-type atomic medium 
the signals described above. The general physical fea- 
ture underlying the mathematical property of complete 
integrability is a delicate balance between the workings 
of dispersion and nonlinearity inherent in the medium. 
Provided that this balance is observed and the system 
is completely integrable it is a general fact that virtu- 
ally any sufficiently intensive localized initial condition 
creates solitons. The overall picture of nonlinear dynam- 
ics can be roughly described as follows. The evolving 
signal created by the incident pulse in the course of the 
dynamics breaks down to a number of solitons and a 
decaying tale. The latter vanishes in due course. The 
soliton-like signals survive (ideally, i.e. in the absence of 
dissipation) for infinitely long time. When the physical 
conditions underlying the complete integrability of the 
optical system are met, the general picture of the nonlin- 
ear dynamics is similar to the described above. Namely, a 
fairly arbitrary localized and intensive initial signal cre- 



ates in the course of nonlinear dynamics a number of 
solitons. The number of solitons can be derived from the 
analysis of the corresponding zero-curvature representa- 
tion. The signal Eq.JTJJ is a rather generic soliton-like 
solution of the nonlinear system Eqs.®, 10 and there- 
fore it is very plausible that such signal can be created. 
Of course an experiment accompanied by a mathemati- 
cal scrutiny would be necessary to further support this 
claim. The purpose of this paper is merely limited to 
pointing to a possibility of interesting physical applica- 
tions. We want to emphasize that in our considerations 
the distinguished role is assigned to the background field 
f2o that turns out to be a nonlinear analog of the con- 
ventional controlling field appearing in the linear EIT 
formulation. The difference between the linear and non- 
linear cases lies in the fact that in the nonlinear case the 
control field and the soliton solution are present in the 
same channel in inseparable fashion of nonlinear super- 
position. In our future work (ArXiv: quant-ph/0411148 
quant-ph/0411149) we plan to further support our results 
on memory reading and transparency control by an inves- 
tigation of transient regimes, when the background field 
Sl adiabatically vanishes or adiabatically approaches the 
point £q. The treatment of relaxation processes in the 
system will be also included. This forthcoming inves- 
tigation will use the method of collective variable |2lJ 
along with other perturbation techniques developed for 
the systems integrable by the inverse scattering method. 



Acknowledgements 

AR wishes to thank for kind hospitality Professor Ulf 
Leonhardt and the University of St. Andrews, Scotland, 
where a part of this work was carried out. IPV acknowl- 
edges the support of the Engineering and Physical Sci- 
ences Research Council. 



7 



[1] L. N. Hau, S. E. Harris, Z. Dutton, and C. H. Behroozi, 

Lett, to Nature 397, 594 (1999). 
[2] C. Liu, Z. Dutton, C. Behroozi, and L. Hau, Lett, to 

Nature 409, 490 (2001). 
[3] D. Phillips, A. Fleischhauer, A. Mair, and R. Walsworth, 

Phys. Rev. Lett. 86, 783 (2001). 
[4] M. Bajcsy, A. Zibrov, and M. Lukin, Lett, to Nature 426, 

638 (2003). 

[5] D. Braje, V. Balic, G. Yin, and S. Harris, Phys. Rev. A 

68, 041801(R) (2003). 
[6] E. Mikhailov, V. Sautenkov, Y. Rostovtsev, and 

G. Welch, J. Opt. Soc. Am. B 21, 425 (2004). 
[7] A. V. Turukhin, V. S. Sudarshanam, M. S. Shahriar, and 

et al., Phys. Rev. Lett. 88, 023602 (2002). 
[8] M. Bigelow, N. Lepeshkin, and R. Boyd, Science 301, 

200 (2003). 

[9] M. Soljacic and J. Joannopoulos, Nature Materials 3, 213 
(2004). 

[10] L. Faddeev and L. Takhtadjan, Hamiltonian Methods in 
the Theory of Solitons (Springer, Berlin, 1987). 

[11] F. Hioe and R. Grobe, Phys. Rev. Lett. 73, 2559 (1994). 

[12] R. Grobe, F. Hioe, and J. Eberly, Phys. Rev. Lett. 73, 
3183 (1994). 

[13] Q.-H. Park and H. Shin, Phys. Rev. A 57, 4643 (1998). 
[14] J. Byrne, I. Gabitov, and G. Kovacic, Physica D 186, 69 
(2003). 

[15] A. Rybin, V. Matveev, and M. Salle, Inverse Problems. 
4, 173 (1988). 

[16] A. Rybin, J. Phys. A: Math, and Gen. 24, 5235 (1991). 
[17] A. Rybin and J. Timonen, J. Phys. A: Math. And Gen. 

26, 3869 (1993). 
[18] V. Matveev and M. Salle, Darboux Transformations 

and Solitons in Springer Series in Nonlinear Dynamics 

(Springer Verlag, 1991). 
[19] M. Fleischhauer and M. D. Lukin, Phys. Rev. Lett. 84, 

5094 (2000). 

[20] A. Its, A. Rybin, and M. Salle, Theor. Math. Phys. 74, 
20 (1988). 

[21] R. Scharf and A. Bishop, Phys. Rev. A 43, 6535 (1991). 



V. APPENDIX. DARBOUX-BACKLUND 
TRANSFORMATION FOR THE A-SYSTEM 



In this appendix we describe the Darboux-Backlund 
(DB) transformation for the A-system. It is plain to see 
that the linear system Eqs.©,© is covariant with re- 
spect to the following DB dressing transformation 

H I = H I -±[D,a 1 {Q)], p = a 1 (A)par 1 (A) (23) 
* = - (71(0)*, (Ti(A) = *i (£i - A) tfT/ 1 . (24) 

Here f is a matrix consisting of three linearly inde- 
pendent solutions of the linear system Eqs.JHJ),© cor- 
responding to three (not necessarily different) values of 
the spectral parameter A. The matrix spectral parameter 



C in our case is defined as 



£ = 




where A', A", A'" are certain values of the spectral pa- 
rameter A. The matrix 4"i is a specification of \& corre- 
sponding to the following particular value of the matrix 
spectral parameter: 



AS 
A* 
A 



We denote the fundamental matrix of the linear system 
Eq.© for A = A as $ . 

It can be shown that for the value of the spectral pa- 
rameter A = Xq the fundamental matrix is i>o = ($q 1 )^. 
Since the subspace of solutions corresponding to A£ is two 
dimensional, the matrix \&i = (ipi,^,^) is constructed 
as follows. The vector ip 3 = ci^ 1 - 1 + c 2 $ 2 ^ + c 3 $[, 3) is 
a general solution of the linear problem with A = Ao. 
Here i = 1,2, 3 denotes a column in the matrix $o- 
To satisfy the structure of the operator a\ Ea. (|24|) we 

require that (^3,^1(2)) = and the vectors '01(2) corre- 
spond to A — Xq. Here, (•, •) denotes a scalar product of 
two vectors in 3D complex vector space. It follows from 
the definition of the matrix $ that (<&q\<I>^) = 5 it j. 
Therefore, we can easily find two appropriate orthogonal 
vectors 01,2 : 



-ci($[, 2) +^ 3) ); & 



- r *<T>( 3 ) 
C 2 V ■ 



The algorithm of finding new solutions of the nonlinear 
system Eqs. can be recapitulated as follows. Find 
a solution <& of the associated linear system Eqs. ©,10, 
corresponding to a certain " seed" solution of the nonlin- 
ear system Eqs.JBJ,©. Build "fi, build a\, use then the 
dressing transformation Eq. 12: ill . 

In the case considered in the present work the matrix 
$0 reads 



/ — tan(r;) 



cos(?7) 



$n = 



— Ao+i-\/— Ag— Q% 
On sin(77) 

-Ao+V-Ag-Og 

1 



Up cosfa) \ 

Ao+i^-Ag-ng ' 
fio sin(^) 

Ao+V~ A o-^o 

1 J 



g l 2(A -A) 



where 



/'1 



exp(ni) 

exp(-/x 2 ) J . r_>r>) 
exp(^ 2 ) 



fee 



I ,• vp-AkA f 

A -A ) t ^(Ao-Al 1 " 



M2 = 



An — A 



8 



The constants 01,2,3 above define the spatial position 
and phase of the slow and fast solitons at a fixed moment 
of time. For simplicity, we assume these constants to be 
real. 

To access the solution Ea.H12jl one has to specify 
ci = ai^o ^fw^nrj, c 2 = a 2 y / e - e\ - ^0, and 



C3 = 0.3 y £0 + V e o — ^0) while A = ieo. This con- 
cludes the construction of the Darboux-Backlund dress- 
ing transformation for the nonlinear A-model. 



